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We study cosmological perturbations in general infiation models with multiple scalar fields and 
arbitrary kinetic terms, with special emphasis on the multi-field extension of Dirac-Born-Infeld (DBI) 
inflation. We compute the second-order action governing the dynamics of linear perturbations in the 
most general case. Specializing to DBI, we show that the adiabatic and entropy modes propagate 
with a common effective sound speed and are thus amplified at sound horizon crossing. In the small 
' sound speed limit, we find that the amplitude of the entropy modes is much higher than that of 

, the adiabatic modes. We also derive, in the general case, the third-order action which is useful 

■ for studying primordial non-Gaussianities generated during inflation. In the DBI case, we compute 
' the dominant contributions to non-Gaussianities, which depend on both the adiabatic and entropy 

M ' modes. 
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CN ; I. INTRODUCTION 

\ While inflation has become a standard paradigm with which to describe the physics of the very early universe, the 
"*7^ • nature of the field(s) responsible for inflation remains an open question. The last few years have seen an intensive 
Qh' effort devoted to trying to connect string theory and inflation (for recent reviews, see e.g. [l|, [3, H, 0, @)i with the 
^ . hope that future cosmological observations, in particular of the CMB anisotropics, could detect some specific stringy 
signatures. 

Of particular interest are scenarios based on the motion of a D-brane in a higher-dimensional spacetime. Since the 
dynamics of a D-brane is described by the Dirac-Born-Infeld (DBI) action, characterized by a non-standard kinetic 
^ \ term, inflation can occur with steep potentials, in contrast with usual slow-roll inflation. In this sense, this scenario, 

■ called DBI inflation 0, H, @ , belongs to the more general class of fc-inflation models [l^, characterized by a 
I Lagrangian of the form P{X, (/)), where X = —d^(j)d^4'l'^- 

, In DBI inflation, the effective four-dimensional scalar field corresponds to the radial position of a brane in a higher 
dimensional warped conical geometry. For simplicity, the other possible degrees of freedom of the brane, namely the 
angular coordinates, are usually assumed to be frozen. Relaxing this assumption and allowing the brane to move 
in the angular directions (see e.g [13- [III) leads to a multi-field scenario, since each brane coordinate in the extra 
\ dimensions gives rise to a scalar field from the effective four-dimensional point of view. 

Beyond the multi-field extension of DBI inflation, it is interesting to study, in the spirit of fc-inflation, very general 
multi-field Lagrangians of the form 



O 

oo 



> 
X 



^' P = P{X",^^), (1) 



with 



X'' ^-\d,^'d>^cj,\ (2) 



where I — I, . . . , N labels the multiple flelds. (In the following we will adopt the implicit summation convention on 
both field and space-time indices.) 

A more restrictive class of models consists of Lagrangians that depend on the global kinetic term X = GjjX^'' 
where the functions Gjj{(j)^) are the components of some metric defined in field space [l^. While this simpler class 
of Lagrangians is enough to describe the homogeneous dynamics of multi-field DBI infiationj_Jt turns out that the 
full inhomogeneous dynamics cannot be described by such a Lagrangian, as we pointed out in [17| and show below in 
more detail. 

The purpose of the present work is thus two-fold. Our first aim is to derive the equations governing cosmological 
perturbations in the generalized models of the form given in Eq. ([1]). Our second aim is to apply this general formalism 
to the multi-field extension of the DBI scenario. 

The structure of this paper is the following. In the next section, we first consider the multi-field DBI action which 
motivates our subsequent study of the general formalism. In section IIIII we derive, in the general case, the field 
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equations for the metric and for the scalar fields, after which we specialize to the homogeneous background. Section 
IIVI is devoted to the dynamics of the linear perturbations in the general case: we derive the second-order action and 
analyse the resulting equations of motion for the perturbations. We then focus, in section |Vl on the specific example 
of the DBI action: we show that the adiabatic and entropy modes propagate with the same speed of sound Cs and 
we compute the second-order action for linear perturbations. For two-field DBI inflation we also compute the power 
spectra of the adiabatic and entropy modes. Finally in section IVTl we discuss non-Gaussianities. We first derive, in 
the general case, the third-order action for perturbations. We then limit our analysis to two-field DBI models, for 
which we compute the main contribution to non-Gaussianity in the limit of small Cg. We summarize our main results 
in the last section. 



II. THE MULTI-FIELD DBI ACTION 



In this section we motivate our reasons for looking at Lagrangians of the general form P{X^'^ , cj)^) by showing, in 
particular, that multi-field DBI inflation is described by a Lagrangian of this form. We also discuss the properties of 
the higher order terms in derivatives which appear in the DBI Lagrangian. 

Consider a D3-brane with tension T3 evolving in a 10-dimensional geometry described by the metric 

ds' = h-''\y'')g^,dx^^dx- + h^/\y'')Gu{y'')dy'dy' = HAsdY^dY^ (3) 

with coordinates Y^ ~ {x^ ,y^}, where /i = 0, . . . 3 and / = 1, . . . , 6 (the label / has been chosen in this way as below 
it will label the multiple effective scalar fields). The kinetic part of the DBI action. 



L=-T3^-dctj^, (4) 
depends on the determinant of the induced metric on the 3-brane, 

7m. = HABd^Y^i-,d,Y^^^ (5) 



where the brane embedding is defined by the functions Y^,^{x'^), with the a;^ being the spacetime coordinates on the 
brane. In our case, they coincide with the first four bulk coordinates. On writing Y^,^ = {x'^ , ip^ {x^^)) , we find 

7m- = (5m- + h Gijd^ip'd.ip^) , (6) 

which after substitution into (|4]) implies 



L = -nh-'^ ^det{Si^ + hGijdi^ip^d.^J) . (7) 
Finally, in order to absorb the brane tension T3, it is convenient to rcscalc in the following way: 

/ - I- , 0' = Vn^'- (8) 
As a result, in the following, we consider the DBI Lagrangian 

P = -j^^{VV-l)-Vi^^) (9) 

where 

V - dct(<5(; + / Gijd^^'d,(j)J) , (10) 

and where we have also included potential terms, which arise from the brane's interactions with bulk fields or other 
branes. From now on we let /= 1, ... , N. 

One can express the above Lagrangian in ^ explicitly in terms of the X^'^ defined in by rewriting V, which 
is the determinant of a 4 x 4 matrix, as the determinant of an iV x matrix: 

V = det((5/ - 2fX/) (11) 

where 

A/ = GikX""-". (12) 



3 



Throughout this paper field indices will always be raised and lowered with the 'field metric' G/j = Gij{(j)^). The 
equality between the expressions pop and pT|) for the determinant can be proved by using the identity det(Id + a) = 
exp[Tr(ln(Id + a))] for the matrix a of components = fGijd^(f)^dy4>^. Indeed from Eq. (fTO|) wc have 



V 



exp 



and on noting that 



Tr(a") = Tr[(-2/A')"] 



(13) 



(14) 



where X represents the matrix of components Xj, one obtains the expression given in Eq. (jlip . 

Another very useful expression for 2? can be obtained by computing directly the determinant in Eq. pop . As we 
show in Appendix 1, this yields 



V = l- 2fGijX" + ^fx'^lxf - 8fx\^XjX^^ + 16f*x\'xjxj^X^^ 



L ' 



(15) 



where the brackets denote antisymmetrisation on the field indices. Wc note that, in four spacctimc dimensions, 
Eq. (jl5p is automatically truncated at order even if the number of scalar fields is larger than four (see Appendix 
1). To use shorter notations, we will rewrite the above equation as 



with 



and where !F{X 



P = 1 - 2fX, 

X = X + J^{X",(b^), 
X = GijX'-' 



IJ 



) can be read from Eq. (fT5l 



T{X 



= -2fxfxf + Afx\^XjX^^ - SfX^^XjX^Xj;^ 



K] 



vJ vK v^l 



(16) 



(17) 
(18) 



(19) 



For a single field, / = 1, it is straightforward to sec that J- vanishes so that the determinant takes the familiar 
form V = 1 + fd^4>d^(j) (for Gn = 1). Similarly, for a multi- field homogeneous configuration in which the scalar fields 
depend only on time and X^'^ = (f)^ (j)'^ /2, one again finds T = Q because of the antisymmetrisation on field indices in 
Eq. p^ . Thus in this case the determinant V reduces to 



V = l-fGij^ 



(20) 



(In the following a bar denotes homogeneous background quantities.) 

However, for multiple inhomogeneous scalar fields, the terms in T , which are higher order in gradients and have 
not been considered in previous works, do not vanish: we will show later in this paper that they drastically change 
the behaviour of perturbations. Furthermore, since they vanish in the homogeneous background, wc expect them 
to modify only the terms in the perturbation equations which contain spatial derivatives. From this discussion we 
therefore see explicitly that the multi-field DBI action docs not depend only on X = GijX^'' (as has been assumed 
in recent works on multi-field DBI inflation [H, [ll]), but requires a general description of the form P = P{X^'-' , 4>^). 

After this digression on the specific form of the multi-field DBI Lagrangian, in the following section we return to 
the general Lagrangian given in Eq. ([1]). 



III. FIELD EQUATIONS 



We begin this section by deriving the equations of motion for the general action 

5=1/ A^x^ + 2P{X'\ 0^) 



(21) 



where we have set SttG = 1. The energy-momentum tensor can be obtained by varying P with respect to the metric, 
and is given by 



(22) 
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where we have defined 

_ 1 f dP dP \ 
P<iJ> ^ 2 \dX^ ^dX^)= ^<"> ■ ^^^^ 

We use this symmetrized derivative of the Lagrangian P with respect to X^^^ for the following reason: since X^"^ is 
symmetric in / and J, the explicit dependence of P on say X^^ and X'^^ can vary depending on the chosen convention 
and the above definition avoids any ambiguity. The same notation will apply to the derivative of any arbitrary function 
which depends on X^"^ . 

The equations of motion for the scalar fields follow from the variation of the action in Eq. pT|) with respect to each 
scalar field. One finds 

-^94V^p</j>a^0'^) + P/ = o, (24) 

where P j denotes a partial derivative with respect to (j>^ . 

Now we consider a spatially flat FLRW (Friedmann-Lemaitrc-Robcrtson- Walker) geometry with metric 

ds^ = -dt"^ + a{tfd^^ , (25) 

where t is cosmic time. Then the kinetic terms defined in Eq. ^ reduce to 

X" = l^'^' , (26) 



where a dot denotes a derivative with respect to t. From Eq. (|22[) . the energy density can be expressed as 

p = 2P<jjyX''' - P (27) 
while the pressure is simply P, and the Friedmann equations are given by 

H^^^ {2P<i.,yX'' - P) , -X'-'P^ij^ . (28) 

The equations of motion for the scalar fields Eq. ([M)) reduce to 

a^'^ (a'^</J></'-') =^,/. (29) 

On calculating the time derivative and taking into accoimt the terms in (j)'^ contained in ^P<ijy^ the above equation 
can be rewritten as 



P<iJ> + P<iL>.<jK>cl}''(b'' )r + [ 3HP<ij> + P<ij>,Kcl)'' ]r -Pi^O (30) 



where, in analogy with Eq. (j23p . we have defined 



_ 1 / dP<ij> dP<u> \ p ^ 

P<IJ>,<KL> = 2 1 Qj^KL + QXLK ] - ^<KL>,<IJ>- (31) 



Finally we end this section by noting that were the Lagrangian to depend on the X^"^ only through X = GijX^"^ , 
then one would define P(X,(j)^) = P{X^^ ,(f)^) so that P^ijy = P.xGjj. In that case all the above expressions 
would reduce to those of \\m. 



IV. LINEAR PERTURBATIONS IN THE GENERAL CASE 



In this section, we derive the second-order action governing the dynamics of the linear perturbations for the general 
action given in Eq. ((2T|) . As in [2^ , we use the ADM approach, and our calculations are very similar to those of [l^ , 
except that we work with the quantities X^'^ rather than X . 

Starting from the metric in the ADM form 

ds^ = -N'^dt^ + h,j [dx' + N'dt){dx^ + NHt) (32) 



where N is the lapse and TV' the shift, the full action becomes 

1 



N + —(-E'^ + E,,E'^) + 2NP 
N 



(33) 



where ^^^R is the scalar curvature of the spatial metric hij with h its determinant, and the symmetric tensor Eij, 
defined by 



(34) 



is proportional to the extrinsic curvature of the spatial slices (Di denotes the spatial covariant derivative associated 
with the spatial metric hi 



The function P = P{X , (j)^^) in Eq. (|33p depends on the kinetic quantities X^"', which can be decomposed as 



with 



(35) 



(36) 



We now work in the flat gauge so that spatial sections are flat, hi 



'^{t)dij and '^'i? = 0. The Hamiftonian and 



momentum constraints, which follow from the variation of (|33p with respect to the lapse and the shift are, respectively, 

2{N^P - P<i.jWv^) +E^ - E,jE^3 = , (37) 



Di 



-A 



- E5i 



j^p<i.,Wd,ci>' . 



(38) 



To zeroth order (background) , the Hamiltonian constraint simply gives the first Fricdmann equation in Eq. (|28p , while 
the momentum constraint vanishes identically. 

At first order in the scalar perturbations, we write 



N=l + 5N, N, = d,i:, (j)^ = ^\t) + Q\t,x) 



(39) 



where, from now on and when there is no ambiguity, we promptly drop the bars on all the unperturbed fields. Note 
that -0 is related to the standard Bardeen potential 5* by 



(40) 



At linear order, the momentum constraint 



gives 



SN = -P,u> 



The Hamiltonian constraint is more complicated, but a straightforward calculation yields 

9V 



2H[^]= 2A5N + Bij(t>-'Q^ + C/Q^ 



with 



A = P<ij>X" - P - 2X'-'X'^''P. 



1.1 yKL I 



<IJ>,<KL>, 



Bij = P<ij> + 2X^'" P<:i.jy^<:KL> , 



Ci 



-Pj + 2P<KL>,lX 



KL 



(41) 



(42) 



(43) 



Actually, this explicit expression for ip is not necessary in order to derive the second-order action, as the terms involving 
Tp (coming from the matter and gravitational parts of the action) cancel each other. The scalar field perturbations are 
related to a useful geometrical quantity, namely the comoving curvature perturbation TZ (see e.g. 21 1 and (2^'). On 
using the standard definition of TZ which combines the metric perturbations with the perturbations of the momentum 
density for the Lagrangian P = P{X^'\ (j)^), one obtains 
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After these preliminary steps, one can now expand the action ([33]) up to second order in the linear perturbations 
SN, ijj and . As mentioned earlier, the terms involving ip cancel each other. On reexpressing 5N in terms of the 
using the constraint (j4T|) . one obtains, after a long but straightforward calculation, 



^(2) = ^ / dtd^xa^ 



(45) 



where the mass matrix is 



MkL = -PkL + 3X^'''P<NK>P<ML> + ^P<NL>4>'' [2P<ij>^kX" - Pk] 



-JpX'^'^ P^NK>P<ML> \X''^ P^ijy + 2P</j>^<AB>-'^^"^-'i''^^] 



1 



and the mixing matrix is 



^KI — 4>'^ P<IJ>,K — J^P<LK>P<MJ>,<NI>X^'^ X^^'' . (47) 

On denoting the coefficient of the kinetic parts in Eq. (|45p by 

Kij = P<,j> + 2P<MJ>,</K>X^^^, (48) 
we find that the equations of motion for the (in Fourier space) are 

KijQ' + ^P<ij>Q-' + (kij + 'iHKij + n.,! - Uij) Q' + (flKi + MiK + iH^Ki) = . (49) 



The propagation velocities can be deduced from the structure of the first two terms in the above equation. On assuming 
that Kij is invertible, the sound speeds correspond to the eigenvalues of the matrix of components (K^^)^^ P^lj> 
(recall that these arc background quantities). 

For a single scalar field, P<:ij> reduces to P,x, and it is easy to see that the kinetic coefficient in Eq. (j48[) is simply 
K = Pjc + 'i.XPxx which can be identified with p,x according to the relation (P7)) . Hence one recovers the familiar 
result [llj that the effective speed of sound is given by 

cl = ^ = — (single scalar field) . (50) 

P,x P,x + 2XP^xx 

For multiple fields and in the particular case where the Lagrangian is a function oi X ^ GijX^'\ i.e. P = P{X, (f>^), 
it has been shown in [l^ (see also [3, [l^) that the propagation matrix (K^^Y^ P^Ljy becomes anisotropic: the 
perturbations along the field space trajectory propagate with an effective speed of sound c^, defined as in the single 
field case above ([50)) . whereas the perturbations orthogonal to the background trajectory propagate at the speed of 
light. In the next section, we examine what happens in multi-field DBI inflation for which P = P{X^'^ , (p^). 

V. LINEAR PERTURBATIONS IN DBI INFLATION 

We now focus on linear perturbations in the specific case of multi-field DBI inflation for which the Lagrangian was 
derived in section 2: 

P{X", </,^) = --^ (VP - l) - V{^'), (51) 

where 2? is given in Eq. (jlip or Eq. (fT5|) . According to Eq. (fT6|) . this Lagrangian can also be seen as a function of X, 
and it will be convenient later to use P deflncd by 



1 



P{X", 0^) ^ P{X, </>^) = --^ ( ^1 - 2f{^')X - 1 ) - y(0^). (52) 
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A. Propagation speed 



Before considering the fuU dynamics of the hnear perturbations, it is instructive to concentrate on their propagation 
speed. According to the general analysis given in the previous section, we simply need to calculate -P</j> as well as 
the matrix Ku, defined in Eq. P5|) . 

To do so, it is convenient to use the form of the determinant given in Eq. (jlip . namely 

I? = det(M), M/ = S/ -2fGiKX^'' . (53) 
In the homogeneous background, the components of the matrix M/ reduce to 

M/ - 6/ - 2fXeie\ (54) 
when expressed in terms oi X = GijX^'^ and of the unit vector in field space 

e^EE^, ei^Gue'. (55) 

We will shortly need the inverse of the background matrix M, denoted by G. Its components are given by 

where, in the second equality, we have introduced the projector orthogonal to the vector e^, 

_L/= (5/ - e/e"'. (57) 

Let us also define 

Cs = v/1 - 2/X = (58) 

which, as we show below, is the propagation speed for all perturbations. Note that this definition coincides with that 
given in Eq. ((50)) for a single scalar field (on replacing P{X, 0^) by P{X , cf)^) given in Eq. (|52p ). and that G given in 
Eq. can be rewritten as 

G/ =±/ +^e/e'^. (59) 

Let us now compute P<ij>- The identity S det{M) = det{M)M^^SM implies, using Eq. (|58p . that 

T^<ij> = -2fciGij, (60) 
where G/j = G^Gkj- It then follows from Eq. ([5T|) that 

P<IJ> = -T^V^IJy - CsGjJ, (61) 

where all quantities are evaluated on the background. For the matrix i^/j, one needs the second derivative of P with 
respect with X^'^ . On using the second derivative of the determinant P, which can be deduced from Appendix 1, it 
is straightforward to obtain 



P<IK>,<.JL> — fcs ^GilGk.! + GijGkl — GikGjlJ ■ (62) 
By noting that X^'^Gkl = X/c^, the above equation together with Eq. (pT|) leads to 



1 ^ 1 



Kij = P<ij> + 2P^iK>.<jL>X'^'^ = -Gu = -^P<i.i> ■ (63) 

Hence we obtain the remarkable result that the propagation matrix is proportional to the identity matrix and that 
all perturbations propagate at the same speed, namely the effective sound speed defined in Eq. (j58p . 

Intuitively one can understand this result as follows. Let us return to the DBI action in terms of the embedding of 
a brane in a higher dimensional spacetime, as discussed in Section [III The perturbations we have considered above 
can be seen as fluctuations of the brane position in the higher dimensional background. Since the brane action is the 
world-sheet volume, its fluctuations propagate at the speed of light from the higher dimensional point of view. From 
a 4-dimcnsional point of view, this translates into the speed of sound as a consequence of time-dilation between 
the bulk time coordinate and the brane proper-time (the Lorentz factor is l/c^). 
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B. Second-order action for the perturbations 

We now turn to the full dynamics of linear perturbations in multi-field DBI inflation. In the second-order order 
action (|45p . the mass and mixing terms could be determined by explicit substitution of P given in Eq. (jSip . Here, 
however, wc follow a more direct route and extend the results of (which were obtained for Lagrangians depending 
only on X). To do so, we use P{X, cj)^) defined in Eq. ([5^ and simply identify the new terms which appear because 
the DBI Lagrangian depends on X rather than X . 

Computation of the first and second order variations of X gives (see Appendix 2) 

5^^'^X = (5(2)1 = (5(2)x + fx Uj h'^diQ^djQ-K (64) 

With respect to the second order action of the extra term in Eq. ([64|) only modifies the spatial gradient term, 
while the rest of that action is unchanged. Hence, as in [l6j . we can rewrite the action in terms of covariant derivatives 
Vi defined with respect to the the field space metric Gjj. This gives 



'(2) 



(65) 



where we have substituted = l/cg and Pxj = fjX/<^l into the expression of [T^, and introduced the time 

covariant derivative VtQ^ = + ^\k^'^ where F is the Christoffel symbol constructed from G/j (and TZiklj 
will denote the corresponding Riemann tensor). Finally the mass matrix which appears above, and which differs from 
Mij in Eq. is 



Mi J = -ViVjP- P^UiKLj 



IK IL 



XP 



H 



■''{P,xj4>i + P.xAj)+ ■ 



1 



P - IH 

2iJ - ' + - 



,2\2 



H 



2iJ2c4 



2c, 



4c3 



2Hci 



(1 + c2) ^i4>j 



2cif^H 



■f,{i(t>j) 



where in the second equality we have substituted the explicit DBI Lagrangian, and used H 
c2 = 1 - 2fX. 



(66) 

-X/cs as well as 



C. Two-field DBI 

We can gain a better intuition for the system of perturbations described by the action by restricting our 
attention to a two-field system, / = 1,2. Then one can unambiguously decompose perturbations into (instantaneous) 
adiabatic and entropic modes by projecting respectively, parallel and perpendicular to the background trajectory in 
field space. In other words, we introduce the basis {co-, £«} where = e^, and is the entropy unit vector orthogonal 
to ei- 



el ^ e\ 



Gije[ei = 1, 



Gijeiei - 0. 



(67) 



We also define 



'2X. 



(68) 



One can reformulate the background equations of motion Eq. (PO]) for DBI in terms of these quantities. The adiabatic 
component is 



CsP.a ~ Cs&^P^ 

whereas the entropy component gives the time variation of e^: 



Xa 



CsP.s I 

— —e, . 



(69) 



(70) 
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In the above equations, P is given in (|52p and partial derivatives with respect to cr or s denote the projection of the 
field space gradients along or el respectively. For example, P^- = e^Pj and P-ss = ^i^i'Di'DjP. 
On introducing the decomposition 

Q' = Qaei + Qsei , (71) 
the equations of motion for and Qs follow from Eqs. ((65)l and ([66]) (see [l^l)- For the adiabatic part one finds 



a4 Qa = (SQ,)- 



[Heir 



(72) 



where the coupling ^ between the adiabatic and entropy components is 



2i 

a 



(l + cj)P,-cfP^,a 



•2 



1 -c? 



P 



/,, + (l+c2)Ks 



while the effective mass of the adiabatic modes can be written in the form 



Ma 



3H -3 



, i&/Hy\ (a/H) 



<j/H J &/H 



The equation of motion for the entropy part can be expressed as 



3H 



Qs 



"2 

2 + + ^ 



Qs 



H 



(73) 



(74) 



(75) 



where the right hand side depends on the Bardeen potential ^P, introduced in Eq. (j40p . and which depends on Q^ and 
Qs through Eqs. (|42)) and ((4T|) . The effective mass appearing above is given by 



1 P^ . . 

-CsP,ss + ^ar^na - + '^clPj.^P,, 



4(l + c,)/3-'^^ 



(2 + c,)(l-c,) 



(1-Cs)' 



1 



l-cV'' Ml + Cs)P^^^ (1 + c,)/ 2/2 ^;-+2'''^^ 



(76) 



{TZg is the scalar Riemann curvature in field space.) Note that in this form, Eq. ((75|) is useful on large scales when 
the right hand side can be neglected — in this case one sees immediately that the entropy perturbation Q^ evolves 
independently of the adiabatic mode. 

In order to study the quantum fluctuations of the system, it is convenient, after going to conformal time r = 
J dt/a{t), to work in terms of canonically normalized fields given by 



3/2 
Cs 



Qa 



(77) 



Remarkably, in terms of these new variables, the second order action (j65p reduces to the very simple form 

1 



'(2) 



I / dr d^T{<2 + „;2 _ 2et;>, - d [{dv^)' + {dv^r 



(78) 



where a prime denotes a derivative with respect to conformal time. The coupling between Vg- and Vs depends on 

e = f 5 (79) 
and we have introduced the two background-dependent functions 

(80) 



aa 



He] 



3/2' 
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This result is similar to that of [1^1 , except for the spatial gradient terms which have the same coefficient for both 
the adiabatic and isocurvature perturbations. The equations of motion for v„ and Vg are 

< - iv'g + [elk' ^-^vs = , (81) 

< + C< + (c]k^ - — + a^A Vs - -iv, = . (82) 
\ a J z 

In the following we will assume that the time evolution of H , a and Cs is very slow with respect to that of the scale 
factor, so that z" / z ~ a" /a ~ 2/t'^. Since r varies from — oo to 0, the wavelength of a given mode is first inside 
the sound horizon (when Ifcc^rl ^ 1) and then crosses out the sound horizon. As in standard inflation, the initial 
conditions for the perturbations are determined by choosing the familiar Minkowski-like vacuum on very small scales. 
Below, we consider in turn the quantization on sub-horizon scales and then the classical evolution on large scales. 



1. Quantization 

For simplicity, we assume that the coupling ^ is very small when the scales of interest cross out the sound horizon, in 
which case one can quantize the two degrees of freedom independently and solve analytically the system (otherwise, one 
can resort to numerical integration by starting deep enough inside the sound horizon as in [H, [^J]). The amplification 
of the vacuum fluctuations at horizon crossing is possible only for very light degrees of freedom. Consequently, if 
/i^ is larger than , this amplification is suppressed and there is no production of entropy modes. Interestingly we 
see from Eq. (|76p that the term coming from the second derivative of the potential along the entropy direction is 
multiplied by the sound speed c^, which implies that, with a similar potential, it is easier to generate entropy modes 
in DBI inflation than in standard slow- roll inflation. Moreover, the second and third terms in ^re always negative 
and thus tend to destabilize the entropic direction. Below we assume that \y,1\/H'^ <^ 1. 

Following the standard procedure (see e.g. [2l| or one selects the positive frequency solutions of Eqs. ([5T|) and 
([82)) . which correspond to the usual vacuum on very small scales: 

v^k^v^k:^ (l - ^) . (83) 

As a consequence, the power spectra for and Vs after sound horizon crossing have the same amplitude 

P..=P. = ^K.P-^. (84) 
However, in terms of the initial fields Qa and Qs, one finds, using ((77)) , 

(the subscript * indicates that the corresponding quantity is evaluated at sound horizon crossing kcs ~ aH) which 
shows that, for small Cs, the entropic modes are amplified with respect to the adiabatic modes: 

Qs* ^ —■ (86) 

Cs 

In order to confront the predictions of inflationary models to cosmological observations, it is useful to rewrite the 
scalar field perturbations in terms of geometrical quantities, such as the comoving curvature perturbation. The latter 
is related to the adiabatic perturbation by the expression (|44)) . which yields 

n = ^Q, , (87) 
a 

so that one recovers the usual single-field result [Tlj that the power spectrum for TZ at sound horizon crossing is given 
by 



27r2 22 4,^20-2 87r2ecs' 
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where e = —H/H^. 

It is then convenient to define an entropy perturbation, which we denote 5, such that its power spectrum at sound 
horizon crossing is the same as that of the curvature perturbation: 

S = Cs^Q,. (89) 

(7 

We thus have 

Vtz,=Vs,=V.. (90) 

We stress that our convention for the definition of S is purely for convenience. 

In contrast with the scalar perturbations, the tensor modes are, as usual, amplified at Hubble radius crossing. The 
amplitude of their power spectrum, given by 

Tr ^ (H!) , ,91, 

\ / k=aH 

is much smaller than the curvature amplitude in the small Cs limit. 



2. Evolution on large scales 

In order to determine the observational consequences of single-field infiation models, it is usually sufficient to evaluate 
the amplitude of the comoving curvature perturbation just after horizon crossing. The reason is that the comoving 
curvature perturbation is conserved on large scales for adiabatic perturbations, as is also the case for the curvature 
perturbation on uniform energy density hypersurfaces Ci which coincides with —TZ on large scales. This property 
is simply a consequence of the conservation of the energy-momentum tensor [2^ (this is also true for non-linear 
perturbations [H, |23, ill ) . 

In contrast with the single-field case, the curvature perturbation generally evolves in time, even on large scales, in a 
multi-field scenario [1^ (see also [l^l for a recent analysis with non-standard kinetic terms}. This can be interpreted 
as due to a transfer between the adiabatic and entropic modes, governed by the relation [l6[: 

n = -S+^^^. (92) 
Cs H 

Note that, whereas this relation might be useful during inflation, it is not always relevant, for example at the end 
of inflation and during reheating where a may temporarily vanish. The importance of the transfer depends on the 
specific model under consideration and can be computed analytically only in some simple cases. 
On large scales the curvature-entropy evolution can be approximated by two equations of the form 

7^ « aHS , S « f3HS, (93) 

where in the latter, we have neglected the second-order time derivative in Eq. (|75p . In our case, the coefficients a and 
(3 are given by 



ST] 1 



c^iJ ' 2 2 

where we have introduced the slow-varying parameters 



and kept only the leading order terms in the expression for /?. 

The system of equations ()93|) can be formally integrated (see [s^l) to yield 

7^ A / 1 TTzs\fn 



(96) 



with 



Tss{t*,t) - cxp ( ^ l3{t')H{t')dt'^ , TnsiU.t) = ^ a{t')Tss{U,t')H{t')dt' . (97) 
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Hence the (time-dependent) power spectra for the curvature perturbation, the entropy perturbation and the correlation 
between the two can be formally expressed as 

Vn^{l + Tls)V*. Vs^TlsV., Cns = {nS) ^ TnsTssV. , (98) 

(recall that TZ and S are implicitly assumed to be uncorrelated at sound horizon crossing). 

An interesting question, which depends on the details of reheating and thus goes beyond the scope of the present 
work, is whether the entropy perturbation during inflation can be transfered to some entropy perturbations after 
inflation, i.e. in the radiation phase. If this is the case, then the primordial entropy fluctuations could be directly 
observable, with the interesting possibility that there could be a correlation between the adiabatic and entropy modes 

In any case, one can introduce the correlation angle O, defined by 

sine = (99) 

which can also be seen as a transfer angle, since 

sine = , . (100) 

If e = there is no transfer [Tus = 0), whereas if |e| = 7r/2 [T-jis 3> 1) the final curvature perturbation is mostly of 
entropic origin. The relationship between the curvature power spectrum at sound horizon crossing and its final value 
is thus 

Vn^ = Vncos'Q. (101) 
This implies, on using the tensor amplitude Eq. (j9ip . that the tensor to scalar ratio is given by 

r = ^ = IGeCsCOS^e. (102) 
Viz 

Interestingly this expression combines the result of fc- inflation [ll| , where the ratio is suppressed by the sound speed 
Cs, and that of standard multi-field inflation [sot . 

From the expression of the curvature power spectrum, one can compute the scalar spectral index in the slow- varying 
approximation. We obtain 

d InT^ dT 
nn = -rr-r = "k. + i?,"^sin(2e)^r^ = nn, - a,sin(2e) - 2(3^sin^e (103) 
ci in rC Oil 

with 



and where we have used 



n-R, — 1 = — 2e* — ?7, — s* , (104) 

H:'^ = -TssP.. H:^^^^-a,-TnsP*- (105) 

The observable spectral index, given in Eq. p03p , not only depends on the values of the various parameters at sound 
horizon crossing, but also on the transfer angle e. 



VI. NON-GAUSSIANITIES 



In the simplest models of inflation, primordial perturbations are characterized by a very small amount of non- 
Gaussianity [20|. However, other models, such as single-fleld DBI inflation, are expected to produce significant 
non-Gaussianity (3^ . If ever detected, primordial non-Gaussianity would be a powerful discriminator between various 
early universe models. In order to study non-Gaussianities, one must analyse the perturbations beyond linear order. 
During inflation, primordial non-Gaussianities can arise from the quantum fluctuations at horizon crossing or, in the 
case of multi-field infiation, from the classical non-linear evolution on large scales (see e.g. (ssl. |33|). 

In this section, we concentrate on the primordial non-Gaussianity originating from the three-point function of the 
scalar field fluctuations, which is the main contribution in single-field DBI inflation. Its calculation requires the third- 
order action in perturbations. Below, we first consider the general case — that is models of the form ([1} — and then 
specialise to DBI. 



13 



A. Third-order action: general case 

We follow the standard approach which has been presented in [s^, [H, [sj , considering successively the third- 
order action from the Einstein-Hilbert term and then from the matter part. A similar calculation of the third-order 
action can be found in [3^, but only for the multi-field Lagrangians of the form P{X, 4>^), where X = 6ijX^'^ . 

The third-order action coming from the gravitational part is the same as in the single field case and is given by the 
expression 



with 



where the relation between 5N and the field perturbations given in Eq. (I4ip can be rewritten as 

5N=NaQ'^ (107) 

with the (field-space) vector 

A/'a = ^P<AB>^''. (108) 

By expanding systematically the matter part of the action up to third order, we finally find (intermediate steps can 
be found in Appendix 3) 

dtd'xa^ {{gi)ABcQ^Q''Q'' + {92)abcQ''Q''Q^ + {9z)abcQ^Q''Q'' 
+ {9a)abcQ^Q''Q'' + {9a)ABQ^d,Q''{5Nr) + {g^)ABQ^d,Q'' {5N^) 
+ {9c)ABcQ^{h''d,Q''d,Q°) + {9d)ABcQ^{h''d,Q''djQ^)] (109) 

{91) ABC = \P.ABC + ]:P.BC.N'a- P<IJ>X^-'NA^fB^fc 
D 2 

+P<I.J>MX'-^AfAAfc - AP<IJ>,<KL>X''^ X'^'^MaMbMc - X'-'P^jjy^BC^A 

+2P<ijy^<KL>,cX'-' X'^'^NaMb 

-\p<IJ>,<KL>,<MN>X"X'^'^X^'''MA^fBMc (110) 

[92) ABC = 2HNaMbMc + \p<KC>,AB(t)^ 

+Ma^'' [5MbX"P<IJ>,<KC> - P<KC>,B + 2MBX'''X'''^P<ijy,<NL>,<KC> 

-2X"P<ijy^<KC>M (111) 

{93)aBC = -Ij^M'APKBO + 2^<BC>,A — ^A [iX^'^ P^IB>,<KC> + X^^ P^BC>,<KL>\ 

-~'2MaP<ij>,<kb>,<mc>X^'^ X^^^ + P<ib>,<kc>,aX^^ (112) 

{9i)ABC = -^4''^' P<BC>,<MA> + -^X^^ (j)^' P^iA>,<KB>,<MC> (113) 
{9a)AB = 2HNaMb + '2NaX''''4)' P<IB>,<KL> - P<IB>,A (114) 
{9b)AB = -P<AB> - 2X^^ P^iB>,<KA> (115) 

{9c)abc = MaX^^P<bc>.<kl> - \P<BC>A - \p<BC>^fA (116) 
{9d)ABC = -]^<l)'^ P<bc>,<ka> ■ (117) 
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In order to get a flavour for the new effects which could arise in multi-field inflation with non-standard kinetic 
terms, it is instructive to compare the above terms with their counterparts in standard multi-field inflation, such as 
studied in [sl]. Substituting the standard matter Lagrangian 

P = GijX'-' - Gu = Si J , (118) 

the above coefficients reduce to 

(51) ABC = ~^Vabc - ]^V^BcJ^A - XMaMbMc, 

{g2)ABc = '2H^^A^fB^fc, {93)abc ^-^-^aGbc, {94)abc = 0, 

{ga)AB = 2HMaMb, {gb)AB = -GaB, {gc)ABC = -^GbC'J^A, {gd)ABC = 0. (119) 

As we will see in the next subsection, the main contribution for DBI inflation will come precisely from the vertices 
associated with the coefficients 34 and gd, which do not exist for standard kinetic terms. 



B. Non-Gaussianities in DBI inflation 



Single-field DBI inflation is an inflationary model which naturally produces a (relatively) high level of non- 
Gaussianity in the small limit, as shown in 0, H^. It is thus important to investigate how the amplitude and 
shape of primordial non-Gaussianities arc modified in the multi-field case [l^, |4l[ . 

Here, we will focus on the dominant contributions to the non-Gaussianities for <C 1, and therefore ignore the 
contributions coming from the gravitational part of the action, which arc known to be sub-dominant. As in the single- 
field case, the dominant contributions come from the terms involving derivatives of P with respect to the X^'^'s, 
because they are enhanced by negative powers of Cs with respect to the other terms. Moreover, terms containing Ma 
are suppressed in the slow-varying approximation. Indeed, if one compares for example the first term of g^ with gd^ 
one finds schematically 



5c ^2 



Hgd H'cs 



(120) 



which is proportional to e = —H/H^ and thus small. Finally, the dominant contributions come from the following 
terms in <S'(3): 

{g4)ijKQ'Q''Q'' + {gd)ijK Q'h^'^d.Q' dkQ'' . (121) 

On substituting the multi-field DBI Lagrangian into the expressions given in Eqs. (|113[) and (I117p . the coefficients 
(74 and gd can be calculated explicitly (one needs the third derivative of P with respect to the X^'^, which can be 
deduced from Appendix 1). 

In the two-field case, decomposing the fields in terms of their adiabatic and entropic components according to (j7ip . 
as well as using (j62p . one finally finds that the relevant terms of the third-order action are given by 

^(main) ^ J ^3 ^ [{Q af + a{Q sf] - ^ [Q. (VQ,)^ - C^Q, (VQ,)^ + 20^4 VQ.VQ, )] |(122) 

where we have used the fact that / ~ l/o'^ in the limit ^ 1. All the terms which appear in Eq. p22p are of the 
same order of magnitude, since Qs — Qa/cs as we have seen earlier. Note that using X instead of X in the DBI 
action (that is, neglecting the higher order terms appearing in J-) would lead to a different third-order action. 

Let us now compute the contribution of these vertices to the relevant three-point functions, by following the 
procedure outlined in detail in [sl]. Working at leading order in the slow varying regime, we use the adiabatic and 
entropic propagators defined by, respectively, 

(Q.(0)g.(r)> - ^(1 - *A:c.r)e'^='^ , (g.(0)Q,(r)) = ^(1 - *fcc,r)e'^-^=^ , (123) 

which correspond to the Fourier transforms of the Green functions, solutions of Eqs. ([5T1) and (|82p with ^ = and 
z" / z = a" /a = 2/r^. The calculation of the three-point functions involve time integrations and we assume that, as 
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usual, the main contribution to these integrals comes from the period around horizon crossing |42l |. which enables us 
to extrapolate the integration bound to t = 0. We also implicitly ignore the correlations at different times between 
the adiabatic and entropy modes, since these are expected to be small if the coupling ^ is small. The quantities a and 
Cs will be considered as constant in time in the integrals. Given these assumptions, the only integrals required are 



iKc 



dr ' 



1 



dr r e 



2AKcs 



2i 



which have be computed by using the appropriate contour in the complex plane (r — > —[oo — ie)). 
The contributions to the three-point function {Qa{ki)Qa(k2)Qa{k3)) are, respectively. 



3H* 



1 i.2l,2jL2 

J. rxi-^ rxj^ A/g 



(124) 



(125) 



from the vertex proportional to and 



K 



'1 



perm. 



(126) 



from the vertex proportional to Qa^yQaY , where we have introduced K = ki + k2 + fcs and used a — H^2tCs 
Summing these contributions, one thus finds 



{QAki)Qa{k2)Q.{kz)) 



1 



kjiki ■ fc2)(2fcifc2 - ksK + 2K^) + perm.] 



(127) 



where the 'perm.' indicate two other terms with the same structure as the last term but permutations of indices 1, 2 
and 3). This is the standard result from single field DBI inflation [39| . 

Let us now turn to the new terms which arise from the entropy fluctuations. They appear in the three-point function 
{Q (j{ki)Q s{k2)Q siks)) , with the contribution 



{2nfS{J2k.) 



1 jU2jL2p 

J. r\j-^i%'2 rug 



(128) 



from the vertex proportional to QaQ^i the contribution 



{27:)H{J2k,) 



1 



fcf(fc2-fc3) ( 



2k2h 
X3 



(129) 



from the vertex proportional to Qa-i^QsY ^^'^ finally the contribution 



(2^)3<5(^fc,) 



i/4 



AV2^eciU^kf 



kl{ki ■ fcz) 



1 fci 



2kik2 



+ {k2 ^ ks) 



(130) 



from the vertex proportional to QsSIQs^Qa- 
Summing these three contributions, we find 



(0.(fci)g.(fc2)Q.(fc3)) = -(2^)-^,5(^fc,) f-,^ /g.-g [-2klklkl 



-kj{k2 ■ fc3)(2fc2fc3 - kiK + 2A'2) 
+kl{ki ■ k2){2kik2 - ksK + 2K^) 
+kl{ki ■ fc3)(2fcifc3 - k2K + 2K'^)\ 



(131) 



As we will see below, the three-point function of the curvature perturbation depends on the symmetrized (with respect 
to permutations of the three wave vectors fc, ) version of this three-point function, and this has exactly the same shape 
as (jl27p . Nevertheless, its amplitude is enhanced with respect to the purely adiabiatic one by a factor of l/c^. 
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Let us now relate the correlation function of the scalar fields to the three-point function of the curvature perturbation 
TZ which is the observable quantity. In order to do so, we use Eqs. ([57)) . ([55]) and ([M]) to write 

with 

X = (f) A = . (133) 

Let us compute the three-point funtion for three wave vectors of comparable magnitude (so that the coefficients 
and As, which depend on the time at which the relevant scales cross the sound horizon, have approximately the same 
value). The two three-point functions of the fields we have calculated give the following contribution 

(7^(fcl)7^(fc2)7^(fc3))''' = {Aaf{QAki)QAi^2)Qaik3)) 

+Aa{Asf{{Qa{ki)Qs{k2)Qs{k3)) + perm.) 
= {A^f{Qa{ki)Qa{k2)Qa{k3)){l + T^s) (134) 

where the adiabatic three-point function is given in Eq. (|127p . Note that the enhancement of the mixed correlation 
{Q a-Q sQ s) by a factor of 1/c^ is compensated by the ratio between Aa and As so that the purely adiabatic and mixed 
contributions in (jl34p are exactly of the same order. 

The superscript (3) in the above equation indicates that we take into account only the contribution from the 
three-point function of the scalar fields. One could also include the contribution from the four-point function of the 
scalar fields, which can be expressed in terms of the power spectra using Wick's theorem, and also from other higher- 
order terms. This has been done for instance in [H, [33|. In the single-field DBI case, the corresponding contibution 
fl^\ is negligible compared to /^]^. Because of the transfer between adiabatic and entropic modes, this should be 
reconsidered in specific multi-field models. Here we simply disregard these contributions though it should be borne 
in mind that they are present in principle. 

Instead of the three-point function, it is now customary to use the non-Gaussianity parameter f^^ defined by 

3 , E^k. 



(7^(fcl)7^(fc2)7^(fc3)) = -(2^)^<5(^ fe,) i^-f^^iv-,^rj ^ . (135) 

From the relation between Vtz and Vtz^, given in Eq. poip . we then obtain, for the equilateral configuration, 

(3) _ 35 1 1 _ _iilcos2e a36) 

One can easily understand this result. The curvature power spectrum is amplified by a factor of (1 + T^^) due 
to the feeding of curvature by entropy modes. Similarly the three-point correlation function for TZ resulting from 
the three-point correlation functions of the adiabatic and entropy modes is enhanced by the same factor (1 -I- T^^). 
However, since /nl is roughly the ratio of the three-point function with respect to the square of the power spectrum, 
one sees that /atl is now reduced by the factor (1 -I- T^g). This may be important in confronting DBI models to 
observations ZsLlil]. 

We end by revisiting the consistency condition relating the non-Gaussianity of the curvature perturbation, the 
tensor to scalar ratio r, and the tensor spectral index nx = — 2e, given in [4^ for single-field DBI. In our case, 

substituting /^]^ ~ — i^cos^O in (|102p . gives 



r + Snr ^ -r i^\J-3 f j , (137) 

As we can can see from (|136l) and (|137p . violation of the standard inflation consistency relation (corresponding to a 
vanishing right-hand side in (jl37p ) would be stronger in multi- field DBI than in single- field DBI, and thus easier to 
detect. In the multi- field case the consistency condition is only an inequality, unless the entropy modes survive after 
inflation in which case O is potentially observable. 
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VII. CONCLUSIONS 

111 this paper we have studied cosmological perturbations in multi-field inflation models for which the Lagrangian 
depends a priori on all the N{N + l)/2 kinetic terms that can be constructed by contracting the spacetime gradients 
of the N scalar fields. Our analysis can be seen as the multi-field extension of fc-inflation, and it also generalizes 
very recent papers which considered more restrictive Lagrangians of the form P = P{X,(j)^). In our very general 
framework, we have computed the second-order action which governs the dynamics of the linear perturbations, and 
were thus able to identify the propagation matrix whose eigenvalues correspond to the generalized propagation speeds. 

We have argued that such a general framework is necessary in order to study multi-field DBI inflation. In that 
model, we showed that all modes propagate with the same speed of sound c^, and hence (if light) they are all amplified 
simultaneously at sound horizon crossing. However, because their respective canonically normalized functions differ, 
the result is that the entropy modes are enhanced with respect to the adiabatic modes: Qs ^ Qa/cs- If there is 
a subsequent transfer from the entropy modes into the curvature perturbation — a generic feature as soon as the 
trajectory in field space is non trivial — the final amplitude of the curvature perturbation is significantly affected by 
the entropy modes. 

We have also derived, in the general case, the third-order action from which one can compute the predictions 
for primordial non-Gaussianities. We have identified the vertices which appear in this action and expressed their 
coefficients in terms of the initial Lagrangian and its derivatives. In the DBI case, we have computed the dominant 
contributions to the non-Gaussianities of the curvature perturbation in the small Cs limit. If there is an entropy- 
curvature transfer, we have shown that the contribution from the entropy modes will increase the amplitude of the 
three-point function with respect to the single field DBI prediction, but the shape of the non-Gaussianities remains 
exactly the same. Since the entropy modes enhance the curvature two-point and three-point functions by the same 
amount, it implies that the /jvl parameter, which is related to the three-point function divided by the square of the 
two-point function, is smaller than in the single field case. The impact of the entropy modes can be expressed simply 
in terms of the entropy-curvature transfer coefficient, which is model-dependent. In the future, it would be interesting 
to study specific scenarios of DBI inflation and to estimate quantitatively this transfer coefficient. 
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Appendix 1: the DBI Lagrangian and its derivatives 
Calculation of the DBI determinant 

The expression for V in Eq. (jlSp can be obtained from its definition in Eq. (jlOp on substituting into 

det{A) = — —ea^a2a3a4.^^^'^'^'^^'^'^ A°'p^A°'p^A°'p^A°'p^ (138) 

the matrix of components 

A''p=S^ + fBfBji, B? = Gijd''(l>', Bl = dp^'. (139) 

On using the identity 

ea,a,030.e"-"^''^>-^^ = -(4 - 3)\j\ . . . Sii (140) 
and the contractions BfB'^ = —2GikX^^, one finally gets the expression 

V = l- 2fGijX^'' + AfX^lxf - 8fx\'x-jX^^ + 16fx\^XjX^X^^ . (141) 

If there arc three scalar fields, the last term disappears because of the antisymmctrization over the field indices. For 
two scalar fields, the last two terms disappear; and for one scalar field, only the first two terms remain. For more than 
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four scalar fields, the truncation at order is natural if one considers T> as the determinant of a 4 x 4 matrix, but 
it is less obvious if one starts from the expression for V as the determinant of a TV x iV matrix, Eq. pTjl . However, 
this can be understood by noting that the term proportional to /" is a sum of products involving n terms of the 
form Xj = If n > 4, among the n terms of the form Bj^, at least two have the same index /.i because the 

spacetime index ^ can only take four different values. Since, by definition of the determinant, all the field indices / 
are antisymmetrized, one thus finds that the term of order /" necessarily vanishes. 

Derivatives 

In order to compute the derivatives of the DBI Lagrangian with respect to X^'\ one can use the explicit expression 
for V given above in Eq. (I14ip . An alternative derivation, which is simpler, is to start from the identity 



V = cxp [Trlnil/] = cxp [Trln(T\f + 5M)] = cxp {Tr [\n{M) + ln(Id + Ar'^SM)] } 

where the components of M, given in (j53[) . arc decomposed into 

M/ - S/ - 2fGjKX^' , SMf = -2fGiK5X^-'. 

Moreover the components of the matrix (M)~^ are the G/ given in Eq. 

Using (|142p . the expansion of 2?^/^ in terms of the matrix U = M~^SM yields 



(142) 
(143) 



^Tr{U) - iTr(C/2) + iTr(C/3) 



1 + iTr(C/) - iTiiU') + i(TrC/)2 + ^TiiU') - ^TiiU)TTiU') + ^{TrUf 



Substituting in the above expression the components 

U/ = -2fGiKSX''' 

of the matrix U, one gets 



(144) 



(145) 



pl/2 ^ 2?l/2 



1 - !Gij8X 



ij 



+f GijGkmGl 



2GilGjk — GijGkl SX 6X 



1 ^ 



N 



-GinGjkGlm — -^GjjGklGmn ) 5X^^ SX^'^^SX^^^ 



(146) 



By interpreting this relation as a Taylor expansion with respect to the variables X^'\ one can obtain the derivatives 
of V^/^ and thus of the DBI Lagrangian with respect to the X^'^ ^ up to third order, as required for the computation 
of the non-Gaussianities. 



Appendix 2: variations of X up to second order 

The computation of the first and second order variations of X follows from 

5^^^X = 5(i)A + .F</j>(5(i)A^'^ + T,kQ^ (147) 

where 

S^^^X'-' = 4>^'Q'^ - 4>'4>-^5N , (148) 

and also from 

-5(2)1 = ,5(2) A + T^u>5^^''X'' + ^^<7,/>,<KL>-5(i)A^'^<5(i) A^^ + I^^klQ^'Q"^ + .F<,,>,k<5(^' A^'^Q^'(149) 
where 

J(2)x^'^ = ^-Q'Q' + \^'^'5N^ + <p^'6^^^v'^ - ^h^'d.Q'd.QJ - 2^^'Q^^6N . (150) 



From the explicit expression for T in Eq. (I19p . we immediately see that its antisymmetric structure implie: 
T^K = P,KL = 0. We also find 

^<ij> = — 2/X _L/j, 

J^<IJ>,<KL> = 2fcl {Gi(kGl)J - GijGkl^ , 

^<ij>j< = — 2/X {Gij^k + Gije^e^^ Glm,k — eie^GjL.K — eje^GiLj<) — '^I.kX _L/j, 
where the first two identities can be deduced from Eqs. (I6ip and (p^ . This readily gives 



Appendix 3: third-order action 

The following relations are useful for determining the third-order action: 



ij 



where we have used Eq. ()108p to rewrite ()148p and (|150|) in the form 



We also have 



P<ij>S'^^^x" = -AP^ijyX'-^ NaMbMcQ^Q^Q^ + GHJVaMbMcQ^Q^Q^ 
'P<bc>MaQ^Q^Q^ + AHNBMAQ^d^Q^SN' 
-P<AB>Q^d,Q^SN' 



as well as 



and 



J(2)p ^ P<,,><5(2)X^'' + \p<u>,<KL>6^'^X'''S^'^X^^^ + P^u>,kQ''5^'^X" + \p,klQ''Q'^ 



S^'^P = P<ij>S^'^X'-' +[p^jj^^^Ki^^S^'^X'-'S^'^X^^ + P^jj^^kQ''S(^^X" 

+ 1 \P<iJ>,<KL>,<MN>5^'^X"S(''>X^^S('^X'^^ + PjjkQ'Q-'Q""' 
6 L J 

+ \ [p<ij>,<kl>,mQ''5^'''X"5'^'^X^^ + P<MN>.ijQ'Q'S^'^X''^ 
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